Gravitational Waves from Quasi- Circular Black Hole Binaries 
in Dynamical Chern-Simons Gravity 
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Dynamical Chern-Simons gravity cannot be strongly constrained with current experiments be- 
cause it reduces to General Relativity in the weak-field limit. This theory, however, introduces 
modifications in the non-linear, dynamical regime, and thus, it could be greatly constrained with 
gravitational waves from the late inspiral of black hole binaries. We complete the first self-consistent 
calculation of such gravitational waves in this theory. We find that future gravitational wave detec- 
tors could place constraints that are 6 orders of magnitude stronger than current ones. 
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Introduction. — General Relativity (GR) has been 
tested to exquisite accuracy in the Solar System and with 
binary pulsars [![, constraining any possible deviations 
when the gravitational field is weak and non-dynamical. 
The non-linear and highly dynamical, strong-field regime 
of GR, however, remains unconstrained. Gravitational 
waves (GWs) emitted during the late inspiral and merger 
of compact objects would be ideal probes of this regime. 
The second-generation of ground-based GW detectors 
(Ad. LIGO 0, Ad. VIRGO Q, KAGRA Q and LIGO- 
India Q) will allow for the first strong-field tests of GR. 

Modified gravity theories that reduce exactly to GR 
in the weak-field, yet deviate in the strong-field, exist. 
One example is dynamical Chern-Simons (CS) gravity @, 
0], where the Einstein- Hilbert action is modified by the 
product of a dynamical scalar field and the Pontryagin 
topological invariant. Such a correction is unavoidable in 
superstring theory to cancel certain anomalies [81-Hoj and 
it also arises naturally in loop quantum gravity, upon the 
scalarization of the Barbero-Immirzi parameter (ill . [l2| , 
and in effective field theories of inflation fl3l | . 

Dynamical CS gravity has not been heavily constrained 
by experiments because (i) it only interacts with grav- 
ity and (ii) it reduces exactly to GR in the weak-field 
limit. This theory has the same post-Newtonian pa- 
rameters as GR in the weak- field [1J, [l5|. The leading- 
order non- vanishing modification to the motion of bod- 
ies enters through corrections to Lense-Thirring preces- 
sion 0. LAGEOS and Gravity Probe B can 
thus constrain the theory, but only extremely weakly. 

GWs emitted during the late inspiral and coalescence 
of black hole (BH) binaries encode invaluable informa- 
tion about the fundamental gravitational interaction in 
the strong field, precisely where CS deviations are largest. 
Supermassive BH mergers are not ideal for such tests be- 
cause their radius of curvature is large, arid thus, CS cor- 
rections are naturally suppressed [19j, |20( ■ On the other 
hand, stellar-mass BHs, the targets of ground-based GW 
detectors, are ideal for testing CS gravity because their 
radius of curvature is small. 



Future GW tests of dynamical CS gravity require the 
calculation of waveform templates with which to filter 
GW data [21| . This is essential because the signal-to- 
noise ratio (SNR) of ground-based detectors is expected 
to be low, at least initially. Therefore, the calculation 
of templates in dynamical CS gravity is urgent, if we 
are to test this theory with future GW observations. In 
calculating such templates, we will discover how a BH 
binary shrinks due to the loss of energy to GW and scalar 
radiation, how the binary's binding energy is modified 
due to the presence of a scalar field and how Kepler's 
third law is corrected. These results are relevant to other 
astrophysical tests of dynamical CS gravity, for example 
with low- mass x-ray binaries pil ] . 

Dynamical Chern-Simons Gravity. — This theory is de- 
fined by the action Q 



S= I d A x\/—g \KgR 



(1) 

where n g = (167rG) , g is the determinant of the met- 
ric g^v, R^vSa and *R^ u P a are the Riemann tensor and 
its dual, R is the Ricci scalar, d is a dynamical field, 
(a, /?) are coupling constants and £ mat is the matter La- 
grangian density. We define the dimensionless parame- 
ter ( = £/m 4 , where £ = a 2 /{K g P) and m is the total 
mass of the system. The characteristic length scale of 
the theory is given by £}/ A and Solar System tests re- 
quire £V4 < C>(10 8 ) km [HI (see for more details). 

Adiabatic Quasi- Circular BH Inspirals. — The inspi- 
ral of comparable-mass compact objects can be studied 
within post-Newtonian (PN) theory, where one assumes 
all characteristic velocities are much smaller than the 
speed of light and gravitational fields are weak [23| . Such 
an approximation at the 3rd post-Newtonian order has 
been shown to be accurate until the very last stages of the 
inspiral (24rl2rj |. We here concentrate on quasi-circular 
orbits, because, by the time GWs emitted in generic or- 
bits enter the sensitivity band of ground-based detectors, 
they will have circularized due to GW emission 27, 28| . 
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This also holds true in dynamical CS gravity [29[, where 
it has also been shown that spinning BHs exist, albeit in 
the slow-rotation limit [3(J HJ . 

A circular orbit is fully described by its binding energy 
E. In dynamical CS gravity, this quantity contains three 
contributions: a gravitational potential energy E v , a ki- 
netic energy E K and a scalar interaction energy E DD . E v 
can be calculated, to leading PN order, via J U[p' 2 d 3 x' , 
where the primes mean that U\ and p 2 are functions of 



Ua is the gravitational potential of BH A [31| 



u A = -2* 

VA 



1 + 3Q 



n A (ij) 



and pa is the density of BH A 

Pa = (m A + ()[[>>,>>,) S^(x k - x k A ) , 



(2) 



(3) 



with the quadrupolc moment 

(201/3584)C(m 4 /mA)xi^i : 
dividual mass, \A 



n 



given by 
311 ] . Here, m,A is the in- 



S\\/mA is the dimensionlcss Kerr 



spin parameter, S\ is the spin angular momentum unit 
vector, and r a and n A are the field point distance and 
unit vector, all relative to the Ath BH, with the angle- 
brackets representing the symmetric and trace-free oper- 
ation, i.e. UA(ij) = nAiTiAj — (^/i)SijnAk'n l A . The poten- 
tial in Eq. ^ can be read directly from the (t, t) com- 
ponent of the metric of an isolated BH in dynamical CS 
gravity 3l|. The density in Eq. ^ must be calculated 
by solving OUa = 4npA, with the potential of Eq. @. 
Combining all these results, 
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where r*i2 and n\ 2 are the binary's orbital separation and 
the separation's unit vector, p = mxm^jrn is the reduced 
mass, and (^4 • B) is the flat-space scalar inner product. 

The scalar field has a rest energy and an interaction 
energy, which is induced because spinning BHs in dy- 
namical CS gravity possess a magnetic-type dipole scalar 
field. When two such BHs are present, the dipole-dipolc 
interaction energy is 



E OD 



25 
256 



(~TXiX2 



(S 1 • S 2 ) - 3(nia • Si)(ni 2 • S 2 ) 



(5) 

This result is derived by analogy with electromagnetically 
charged BHs in GR [32|. The kinetic energy to leading 
PN order is E K = pv 2 /2, with v the relative velocity. 

Let us now re-express E v , E DD and E K in terms of 
u = (7TTO/) 1 / 3 = (tow) 1 / 3 where / is the GW frequency 
and oj is the orbital angular velocity. This is achieved by 
finding the leading-order, effective-one-body, equation of 
motion for the binary constituents, which in the center 



of mass frame is p ri2 oj 2 = F v + F DD . The magnitude of 
the gravitational and dipole-dipolc forces are computed 
by differentiating the potential and dipolc-dipole ener- 
gies with respect to ria- We have checked that in the 
test particle limit, this equation of motion reduces to the 
modified Kepler's third law found in [3l|. From this, we 
obtain Kepler's third law ri2 = (to/u 2 )(1 + 5C r u A ) with 

25 1 / » 

SC r = --C-X1X2 ((Si • S 2 ) - 3(m a ■ Si)(n 12 • S 2 ) 
512 -q \ 
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where we have orbit averaged, as implied by (•••)„. We 
also obtain v = ri 2 u> = u(l + 5G r u A ). 

We now have all the ingredients to compute the binding 
energy in the center of mass frame in terms of u. Com- 
bining E v , E DD and E K and taking the orbital average, 
we find E = -{p/2)u 2 {\ + 8C E u i ) } with 



5C E = ^tC^ ( (Si ■ S 2 ) - 3(m a • Si)(m 2 • S 2 ) 
256 rj 



+ 
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where -q = p/m is the symmetric mass ratio. This re- 
duces to Eq. (100) of [3l| in the test particle limit. 

For a circular orbit, the balance law states that the 
rate of change of the binding energy E must be exactly 
balanced by the flux of scalar and gravitational radia- 
tion taken out to infinity and into any horizons. One 
can then show that E = -{Z2/b)p 2 v i uj 2 [l + (5C B u 4 ] = 
-(32/5)t ? 2 w 10 [1 + (SC E + 4<5CV)u 4 ] where 



SC K 
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(ij)- 



(8) 



Here, v\ 2 = v{ — v 2 is a unit vector pointing in the di- 
rection of the difference of the orbital velocities, v\ 2 = 
v\ 2 v{ 2 , and A* = (m 2 /m)xiS\ — (mi/m)x2S 2 . The first 
term is a CS correction due to the emission of scalar ra- 
diation, while the second one is a CS correction to the 
emission of gravitational radiation. 

GWs from Quasi- Circular BH Inspirals. — The bal- 
ance law allows us to write an evolution equation for the 
GW frequency: / = / G r(1 + SCu 4 ) and to leading order 
in the PN approximation / GR = (96/5)n 8 / 3 M 5/3 f 11/3 + 
0(u 13 ) while SC = SC E - 3SC E + A5C r , or 
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This evolution equation then gives the time-domain rep- 
resentation of the GW phase via 4>(t) = J 2n(f / f)df. 

In the extraction of GWs from noisy data, it is custom- 
ary to employ the Fourier transform of the GW response 
function. This quantity can be computed analytically via 
the stationary phase approximation [33i - i36j . Neglecting 
PN amplitude corrections to the GW response, the sky- 
averaged Fourier transform is h(f) = Af~ 7 / 6 exp[i^(/)]. 
The overall amplitude is the usual GR quantity: A = 
3Q- 1 / 2 ir-' 2 / 3 M 5/6 Dl 1 , where M = ?y 3/5 m is the chirp 
mass and Dl is the luminosity distance. The Fourier 
phase *m = * OR (/) + where tt OR (/) is the GR 

result [37H39| . while 

6*(f) = ^(7rMf)- 5 / 3 {-W6Cu 4 ) (10) 

is the CS correction. This correction enters at 2PN order, 
and thus, it is degenerate with the spin-spin PN correc- 
tion to the GR Fourier phase. 

Such a modified gravitational waveform is naturally 
contained in the parameterized post-Einsteinian (ppE) 
framework [4^], aimed at proposing waveforms that 
model-independently deviate from GR. The simplest ppE 
Fourier phase is * PPE = * GB + /3 PPE (7r.M/) bppE , which 
can be mapped to Eq. (fT0|) through the ppE parameters 



(ii) 



The CS correction depends on the coefficient SC, which 
in turn depends on a complicated expression [Eq. ©]. 
Such an expression simplifies greatly when the spin an- 
gular momentum is aligned or counter-aligned with the 
orbital angular momentum: 



SC 



309845 1 
553728 V 



47953 
61969' 



1 - 



199923 
61969 



(12) 



This result seems divergent in the extreme mass-ratio 
limit SC oc C(Xs — Xa) 2 /v 2 ■ However, if we take 

the limit mi — > 00 while keep ing mi and u fixed, SC 
actually scales as rrC^rn^ 2 [29|. Astrophysical BHs have 
a theoretical lower mass bound of 3Mq, but one might 
want to consider a primordial BH with arbitrarily small 
mass. If we take m 2 — > 0, SC does diverge, but this is 
beyond the small-coupling approximation that we have 
used, as we explain below. 

Validity of Approximations. — We have here employed 
several approximations. First, we assumed a small cou- 
pling, i.e. £ <C M 4 where M is the smallest mass scale 



of the problem. This approximation is required because 
Eq. ([T]) represents an effective theory; an expansion to 
quadratic order in the curvature. Higher curvature terms 
are being neglected and these would enter with higher 
powers of £' = £/Af 4 . Therefore, the field equations de- 
rived from Eq. ([1} are also only valid to linear order in 
C . It is obvious that if one takes one of the BH masses 
to zero, the small coupling approximation is violated. 
Second, an exact, closed- form solution that represents a 
spinning BH in dynamical CS gravity is currently known 
only to 0(xa) m a XA "C 1 expansion. Therefore, the po- 
tential energy in Eq. ([]} is formally only valid to 0{x\)i 
and thus, so is the waveform computed above. 

Since scalar radiation is caused by scalar dipole 
charges, we can estimate the accuracy of the slow- 
rotation expansion by computing the dipole charge to 
all orders in \ and then comparing this to an expression 
truncated to 0{x 2 )- The dipole charge ha is defined as 
the asymptotic coefficient of the r~ 2 cos 9 term in the so- 
lution for when considering an isolated spinning BH in 
dynamical C S g ravity. The variation of the action with 
respect to d |30l l4~fl] leads to the 1? equation of motion, as 
given in Eqs. (10) and (25) of [3l|. Employing a multipo- 
lar decomposition, we can solve this equation at dipole 
(I = 1) harmonic order using Green's function methods 
and exactly obtain the scalar dipole charge as 



a2 



2%i ' 2vT 



Xi-xi(3-2v / l^d) 



2x^ 



(13) 

For \xa\ < 0.8, the difference between /4f lU ' ) anc l ^ s trun- 
cated expansion at 0{x 2 ) is always less than 10%. 

Future Constraints with GW Observations . — Let us as- 
sume that a GW observation has been made and found 
consistent with GR. One can then ask how large £ can be 
to be consistent with such an observation, thus placing a 
constraint on l^ 1 / 4 !- One can estimate such a constraint 
through a Fisher analysis |34| . For sufficiently high SNR, 



the accuracy to which a given parameter 6 a can be mea- 
sured can be estimated via (A9 a ) = ^(T- 1 ) aa (no sum- 
mation implied), where 



r 



ah 



4Rc 



dgh(f)d b h(f) 
Sn(f) 



df-, 



(14) 



is the Fisher matrix, partial derivatives arc with re- 
spect to 6 a , and S n (f) is the noise spectral density for 
Adv. LIGO (also for Adv. VIRGO and KAGRA) [Hj], 
ET 0, LISA [3 and DECIGO/BBO [H. The limits 
of integration / m i n and / max will be here taken to be 

/min = max(/i ow , /i yr ), / max = min(/ high , / end ) , (15) 

where /i ow and /high are lower and higher cutoff frequcn- 

is the 



cies of a given detector, respectively, while /i yr 
GW frequency 1 year prior to coalescence and / enc j is the 
frequency at the innermost stable circular orbit which 
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can be obtained by solving Cq = where Cq is given in 
Eq. (1.5) of 0. 

Let us first concentrate on the spin-aligned and spin 
anti-aligned cases. Such a waveform has eight parame- 
ters: {8 a } = {hiM,lmi,Xs,Xa,t c ,<f>c,DLX}- We here 
choose the fiducial values t c = 4> c — Q = and impose 
the prior \A < 1- ET and DECIGO/BBO should allow 
us to constrain 

£i/4 < o(10 - 100)km. As expected, the 
constraint scales as the smallest length scale of the tar- 
get system (i.e. the horizon size of the smaller BH), and 
hence it is of 0(M). Due to this reason, LISA can only 
place t; 1 / 4 < O(10 5 — 10 6 )km. Second-generation ground- 
based detectors are not sensitive enough to measure ( 
with non-precessing binaries, within the small coupling 
approximation, due to degeneracies between £ and xa- 

Let us now consider spin precessing binaries. In such 
systems, all momenta precess around the total angular 
momentum, thus breaking certain degeneracies between 
the CS corrections and the spin-spin GR couplings. The 
CS correction to the precession equations will also en- 
ter at 2PN order relative to GR, but since the effect of 
spin already enters at 1.5PN order relative to the lead- 
ing Newtonian term, such CS precession corrections can 
be neglected. Therefore, it would suffice to evolve pre- 
cessing spins with GR evolution equations 47|, |48[ . To 
simplify the problem further, instead, we will assume 
that when spins precess, Xi,2 become uncorrelated to the 
other template parameters, which amounts to repeating 
the calculation of the spin-aligned case, but assuming the 
spin parameters are known a priori. This approximation 
has been shown to be very accurate, for example when 
considering GW bounds on the graviton Compton wave- 
length 0]. 

The projected bounds on £ 1//4 for systems with precess- 
ing spins under this approximation are shown in Fig. Q] 
as a function of (m, X2), with (mi/m 2 ,Xi) = (2,0.4), 
and at a fixed luminosity distance of Dl = 0.1 Gpc for 
Adv. LIGO and D L = lGpc for ET, LISA and DE- 
CIGO/BBO. The colored contours show constraints that 
satisfy AC < 1 [where the smallest length scale of the 
system is taken to be the horizon of the smaller BH, i.e. 
M — 7712(1 + y/l — x|)]' as otherwise the small-coupling 
approximation is violated. Of course, the regions that 
satisfy AC < 1 depend on the choice of Dl: the col- 
ored contours would be larger if one detects a GW from 
a closer binary and thus with a higher SNR. We have 
here chosen a rather small Dl (especially for Adv. LIGO) 
in view of the expected detection rate. The bounds one 
could place on t; 1 / 4 are now stronger as precession has de- 
stroyed correlations that are present in the spin-aligned 
or anti-aligned case. This time, there are parameter re- 
gions where second-generation ground-based detectors, 
such as Adv. LIGO, Adv. Virgo, KAGRA, as well as ET, 
could constrain dynamical CS gravity to 



This is six to seven orders of magnitude stronger than 
current Solar System bounds fl6| . 
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FIG. 1: Projected lcr-constraints on C in km with 
Ad. LIGO (top left), ET (top right), DECIGO/BBO (bot- 
tom left) and LISA (bottom right) at a fixed luminosity dis- 
tance (shown on the top right of each panel), spin-aligned BH 
binaries with (Xs,Xa) assumed measured, mi/7712 = 2 and 
Xi = 0.4. This roughly models projected constraints with pre- 
cessing BH binary observations. The colored contours show 
the regions of parameter space where the constraints on £ 1//4 
also satisfy AC < 1, and thus, the small coupling approxima- 
tion is satisfied at the fiducial luminosity distances chosen. 



Future Work. — This paper opens the door to several 
follow-ups. One possibility is to include eccentricity and 
precessing spins in the inspiral evolutions to estimate how 
much these effects improve or deteriorate the projected 
constraints. One could also carry out a Bayesian param- 
eter estimation and a model-selection study to estimate 
projected constraints for signals with low SNR (4^ |. 

The projected constraints presented here rely on sev- 
eral approximations, but all systems analyzed were cho- 
sen such that these constraints are rather conservative. 
Perhaps the dominant assumption is the truncation of 
the inspiral at the ISCO. Indeed, CS corrections become 
larger as the binary shrink, so if one could integrate all 
the way up to the light-ring and merger, the projected 
constraints would be much stronger. For this, one would 
have to develop an effective-one-body resummation of the 
waveform, which would in turn require the numerical sim- 
ulation of BH binaries in dynamical CS gravity. 

One could also investigate highly spinning BHs and 
other neutron star binaries. We have succeeded in ob- 
taining the dipolc charge for arbitrarily rapidly rotating 
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BHs, but the waveform also depends on the deviation of 
the BH metric's quadrupolc moment. A study of binary 
neutron stars in dynamical CS gravity might allow im- 
mediate constraints through binary pulsar observations. 
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